double* dvecalloc(char n);

void dvecfree(double* v);

double** dmatalloc(char nx, char ny);

void dmatfree(double** m, char nx);

void _vcopy3(double* src, double* dst);

void _vcopy4(double* src, double* dst);

double** T3(double** m);

double** mmult3(double** m1, double** m2);

void _mmult3(double** m, double** m1, double** m2);

double** mmult4(double** m1, double** m2);

void _mmult4(double** m, double** m1, double** m2);

double* mvmult3(double** m, double* v);

void _mvmult3(double* p, double** m, double* v);

double* mvmult4(double** m, double* v);

void _mvmult4(double* p, double** m, double* v);

double* vadd3(double* v1, double* v2);

void _vadd3(double* v, double* v1, double* v2);

double* vadd4(double* v1, double* v2);

void _vadd4(double* v, double* v1, double* v2);

void _vscale3(double* v, double s);

void _vscale4(double* v, double s);

void _normalize4(double* v);

void vprint(double* v, char n);

void mprint(double** m, char nx, char ny);
